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^ ; Abstract 
T— I ■ 

^ . Relativistic Hamiltonians are defined as the sum of relativistic one-body ki- 

, netic energy, two- and three-body potentials and their boost corrections. In 

I this work we use the variational Monte Carlo method to study two kinds 

' of relativistic effects in the binding energy of and ^He. The first is due 

OO ■ to the nonlocalities in the relativistic kinetic energy and relativistic one-pion 

. exchange potential (OPEP), and the second is from boost interaction. The 

i OPEP contribution is reduced by ~ 15% by the relativistic nonlocality, which 

i-H I may also have significant effects on pion exchange currents. However, almost 

^ ' all of this reduction is canceled by changes in the kinetic energy and other 

Ch ■ interaction terms, and the total effect of the nonlocalities on the binding en- 

^ . ergy is very small. The boost interactions, on the other hand, give repulsive 

! contributions of ~ 0.4 (1.9) MeV in (^He) and account for ~ 37% of the 

H ' phenomenological part of the three-nucleon interaction needed in the nonrel- 

■ ■ ■ ' ativistic Hamiltonians. 

PACS numbers: 21.45. -hv, 21.60.Ka, 24.10.Jv 
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I. INTRODUCTION 



It is generally accepted that QCD is the fundamental theory of strong interactions, 
however, due to quark confinement, the genuine QCD degrees of freedom are not explicit 
at low energies. In low energy nuclear physics, nucleons and mesons are believed to be the 
physical (effective) degrees of freedom. In the nonrelativistic many-body theory, nuclei are 
regarded as bound states of nucleons interacting via two- and three-body potentials. All the 
sub-nucleonic and meson degrees of freedom, as well as relativistic effects are, in some way, 
absorbed in these potentials. Typically the nonrelativistic Hamiltonian is expressed as 

^^« = E^ + E%+ E ^... + ---, (1-1) 

i '^'"'i i<j i<j<k 

and models of two- and three-body potentials are constructed by fitting observed data. The 
ellipsis in Eq. ( |1 . 1| ) represents A^-body interactions (A^ > 4) which are thought to be much 
smaller than two- or three-nucleon interactions, and therefore neglected. 
The central problem is to solve the many-body Schrodinger equation: 

Hnr\^)=E\^). (1.2) 

The eigenvalues E can be compared with experimental energies, and the eigenstates |\E') 
can be used both to study the nuclear structure and probe it through electron-nucleus 
scattering experiments, and to calculate rates of nuclear reactions which may have important 
applications in several domains of physics. 

Schrodinger equation ( |1.2|) is difficult to solve due to the large spin and isospin depen- 
dence of Vij and Vijk- Several techniques have been developed, among which are Faddeev- 
Yakubovsky Harmonic- Hyperspherical basis 0, and Quantum Monte Carlo (QMC) 
methods. The first two methods are limited to solving 3- and 4-nucleon systems, whereas 
with the third method it is now possible to calculate the ground state energy and wave 
function for A=2-8 nuclei with great accuracy. 

Some of the results obtained by Pudliner et al. are listed in Table I. The Argonne 
^18 2-body potential p, fitted to A^A^ scattering data and the deuteron binding energy, 
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and Urbana IX 3-body potential constrained to give the correct binding energy of 
and density of nuclear matter, are used in these calculations. It works rather well for "^He, 
however, as can be seen from Table I, the A=6, 7 and 8 nuclei appear to be systematically 
underbound. It is interesting to note that a large fraction of the total Vij comes from the 
one-pion exchange potential (OPEP) and the dominant part of Vijk comes from two-pion 
exchange. Also notice that the three-body interaction is much smaller than the two-body 
interaction, yet it is crucial to obtain the observed energies, because of the large cancelation 
between the kinetic energy and the two-body potential energy. 

Although the nonrelativistic QMC techniques have advanced to such a level that the 
binding energies of light nuclei predicted by a realistic Hamiltonian can be calculated with 
<1% error 0J^, the effective description of nuclear dynamics by means of nonrelativistic 
Hamiltonians may have intrinsic deficiencies. In particular, when the nonrelativistic poten- 
tials are fit to the experimental data, relativistic effects are automatically buried in these 
potentials. How well can these effects be represented by means of local nonrelativistic poten- 
tials is an important question to be answered. In other words, we may investigate whether 
an explicit and more correct treatment of relativistic effects can resolve the systematic un- 
derbinding of the nonrelativistic results for A = 6,7,8 nuclei. 

Furthermore, with the recently completed multi GeV electron accelerator facilities such 
as TJNAF, experiments will be performed at energy and momentum transfer regimes where 
relativistic effects are substantial. Clearly, investigation of these effects has become in- 
creasingly important. Above all, no matter how small the relativistic effects might be, 
understanding them is a fundamental quest, just like understanding the fine and hyperfine 
structures in the hydrogen atom. 

Several approaches have been developed to study various aspects of the relativistic effects 
in few-body nuclear physics. They can be classified in two main categories: effective field 
theories and relativistic Hamiltonian dynamics. Within the first one, the Bethe-Salpeter 
equations for the two- and three-body systems have been solved using a separable kernel [0 . 
Also covariant three-dimensional reductions of the relativistic integral equations have been 
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applied, along with one-boson exchange models for the kernels, to the three-nucleon sys- 
tem. Here we refer to minimal relativity in the Blankenbecler-Sugar equations and, more 
recently, to the spectator (Gross) equations 0. In the relativistic Hamiltonian dynamics 
approach, relativistic covariance is achieved through the Poincare group theory. One appli- 
cation of this method is the light front dynamics, which has been applied to the two-body 
system |]10|, and the other is the instant form. 



The present interest in the relativistic Hamiltonian dynamics in the instant form stems 
from the fact that the ground states of the Hamiltonian can be studied with the Quantum 
Monte Carlo methods that have already been developed in the nonrelativistic approach. 
Earlier |1ll| , p!2| and the present study are limited to the A=2, 3 and 4 nuclei, but attempts 



to study larger nuclei are in progress. In addition, the use of short-range phenomenological 
terms in the interaction gives the flexibility to allow a very good fit to the two-body scattering 
data with ~ 1. And finally, it is not obvious that the entire short and intermediate range 
2-nucleon interaction can be represented as due to the exchange of a few types of mesons. 
Thus more general ways of studying relativistic effects in nuclei are desirable. 

In this paper we report new results for the binding energies of the A = 3,4 systems, 
using the relativistic Hamiltonian dynamics in the instant form, where for the first time the 
nonlocalities induced by the relativistic effects in the one-pion-exchange potential (OPEP) 
are taken into account. In Sec. II we discuss the relativistic Hamiltonian used in this work. 
In Sec. HI we apply Variational Monte Carlo (VMC) techniques and present results. Finally 
we summarize in Sec. IV. Some of the detailed derivations involved in this work are given 
in the Appendix. 

II. THE RELATIVISTIC HAMILTONIAN 

In relativistic Hamiltonian dynamics in instant form the momentum (P) and angular 
momentum (J) generators are chosen in the conventional way and therefore are independent 
of interaction, while the Hamiltonian (H) and boost (K) generators have interaction terms. 
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Based on the pioneering work of Bakamjian and Thomas [jT3[ and Foldy [jl4[, the relativistic 
Hamiltonian can be expressed as: 

Hr = Y1 U'^f+Pi - "^i) + E + + E \^i3k + 5^*jfc(Pijfc)] + ■ ■ ■ , (2.1) 



i<j i<j<k 

where Vij are two-body potentials in the "rest frame" of particles i and j (i.e. the frame 
in which Pij= Pi+Pj=0). Similarly Vijk is the three-body potential in the frame in which 
Pijk=Pi+Pj+Pk=0. The 6vij(Pij) and SVijkiPijk) are called "boost interactions" and de- 
pend upon the total momentum of the interacting particles. Obviously, 6vij{Pij=0) and 
6Vijk{Pijk=0) vanish. The Hnr contains approximations to the kinetic energy T, Vij and 
Vijk, and totally neglects the boost interactions. In the case of deuteron, we can always go 
to its center of mass (cm.) frame where total momentum Pjj=0. The Vij is adequate to 
describe the deuteron in its "rest frame", however, in A > 2 nuclei the total momentum of 
any pair of nucleons is not necessarily zero in the cm. of the whole nucleus, therefore the 
interaction between the pair can not be correctly described by Vij alone. 

Two kinds of relativistic effects in the interaction are studied in this work: the boost 
interaction Svij{Pij) due to the motion of the cm. of nucleons i and j in the rest frame of 
the whole nucleus, and the nonlocality due to the relative motion of two nucleons in their 
own cm. frame. The latter will affect the "rest frame" potential Vij. The boost interaction 
has been studied in detail in Refs. |[T7|JTT|JT^ , so we will only give a brief discussion. On the 
other hand, the treatment of the nonlocality of OPEP in quantum Monte Carlo calculations 
is discussed in detail. 

The boost interaction 6v{Pij) is determined from the "rest frame" potential "Vij through 
relativistic covariance [|15|,|16[. The 6v{Pij) is expanded in powers of Pfj/Am? and only the 
leading corrections are considered in this work. The 5v(Pij) is given by: 

5t;(P) = --^5 + -^ [P-rP-p,5 ] + -^ [(o-i-^T2)xP.p,C], (2.2) 

where the subscripts ij of v, P, p and r have been suppressed for brevity. Here p=(pi — Pj)/2 
is the relative momentum operator, and cr=2s are the Pauli matrices for spin 1/2 particles. 



Various aspects of (5f(P) are discussed in Ref. [|T^. The first two terms of Eq. ( p.2|) 
are denoted as 5vre and 5fLc; they have simple classical origins in the relativistic energy- 
momentum relation and Lorentz contraction. The last term contains contributions from 
Thomas precession and quantum effects. They are denoted as 5vtp and 5vqm and are much 
smaller than the first two terms. For example, the contributions of 5vre-, ^vlc, Svtp and 
6vqm to the energy of triton are found ^ to be 0.23(2), 0.10(1), 0.016(2) and -0.004(2) 
MeV, respectively. Since the main contribution comes from the first two terms, for simplicity, 
we neglected the last two terms in the 3-, 4-body calculations in this work. 

In addition to the boost interaction, another source of relativistic effects comes from the 
nonlocality. In the following discussion we will use the two-nucleon cm. frame in which boost 
interaction vanishes, and focus on the two-body "rest frame" interaction. In most existing 
nonrelativistic potential models, the OPEP has been calculated using the nonrelativistic 
Pauli spinors. Without ttNN form factors, it is given in momentum space by: 

Vtt,NR[(1) — n 5— o ) I'^-'JJ 

where J^nn is the pion-nucleon coupling constant, fi is the pion mass and q is the momentum 
transfer, 

q = p - p'. (2.4) 

Here p and p' are the initial and final momenta of nucleon i in the center of mass frame, 
and the Vt,^nr is local, i.e. it depends only on q. 

In contrast if one uses relativistic Dirac spinors the on-shell OPEP has the form: 

TTl Tfl 

V "iTi^ + y V ^ + V 

This potential is dependent not only on q but also on p and p', which results in a nonlocal 
potential in the configuration space. The interaction ( p.5| ) is regarded as energy independent 
and used in many-body Schrodinger equations. By expanding the square roots it can be 
easily verified that the leading correction {vn,Rei — ^tt.atr) is of order p'^/m'^, i.e., of order 
f^/c^ where v denotes the velocity of the nucleons in the center of mass frame. 



In Ref. it is shown that the relation between the boost interaction Sv{P) and the 



static vj\ffi is independent of the origin of vnr up to order P'^/m'^, and presumably beyond. 
Thus the knowledge of the static vnr is sufficient to obtain 6v{P). In contrast the above 
relation between Vn,Rei and Vn,NR (Eq- [275| ) is specific for the interaction due to exchange 
of pseudoscalar mesons by Dirac nucleons via either pseudoscalar or pseudovector coupling. 
By expanding the square roots in Eq. (|2.5|) we obtain: 



V.,Rel{P, P) = V7r,NR{c[) " ^ " £^ + " " ' (2-6) 

whereas the interactions generated by exchange of scalar (S) or vector {V) mesons have 
different relations 111 711: 



vsMp'^ p) = vs,nr{^) (l - + ■■■)' ^^-^^ 

vvMp, p) = vv,NRicD (^1 + ^^^^ + ■■■)■ (2-8) 

Realistic models of nuclear forces contain momentum dependent terms which presumably 
take into account some of the relativistic corrections to the phenomenological short and 
intermediate range parts of v^^. However, most configuration space models do not contain 
long range, momentum dependent tensor forces occuring in Vt^^r^i- In exact calculations the 
tensor force can not be generally expanded in powers of /m? . In any truncated expansion 
the force diverges at large values of p, and can yield divergent attraction. 

In the present work the two-nucleon potential is expressed as: 

VnN = V^,Rel + VR (2.9) 

where vr is the remaining part of the 2-body potential which is phenomenological. We can 
also write the OPEP given in Eq. (|2.5| ) as: 

Vn,Rel = V7T,NR + (^Tr.ReZ " Vn,NR) ■ (2.10) 

The term in parenthesis is the relativistic correction. The nonrelativistic potential models 
do not consider this correction explicitly: the data is fit using Vn,NR in Eq. ( ^^ , thus some 



of its effects go into the phenomenological part of the potential vr. The vr in relativistic 
Hamiltonian differs from that in nonrelativistic Hamihonian due to the difference in 'v.,r,Rei 
and Vj^^NR as well as that in the kinetic energy operators. 

We construct our Hr to be phase equivalent to the isoscalar part of the nonrelativistic 
Hnr containing Argonne fig. The relativistic effects can then be studied by comparing 
results obtained from our Hr and the isoscalar Hnr without considering the small isospin 
symmetry breaking terms in the latter. The relativistic Hamiltonian for two-nucleon system 
in its center of mass frame is chosen as: 



Hn = 2v'p2 + rn^-2m + , v^,nr{ci) . , , , + vr. (2.11) 

where vr has the same form as the isoscalar part of Argonne 0: 

~^R= ^p(^^.)Of,. (2.12) 

p=l,14 

Oi~^'^^ = l,Ti ■ Tj,(Ti ■ (Tj, {(Ti ■ (Tj){Ti ■ T j) , Sij, Sij^T j ■ Tj), L ■ S, L ■ S(r j ■ T j) , 

L\ L\t, ■ r,), L\a, ■ a,), L\a, ■ o-,)(t, ■ r,), (L ■ S)^ (L ■ Sf{T, ■ r,). (2.13) 

The 'U,r,Arii used in Argonne fig is given by: 

v.,nr{v) = Ifi^ [Y^ir)(T, ■ + T^ir)S,,] r, ■ t„ (2.14) 
6 47r 



where 













(l-e- 
3 


fir 






— + - 

/ir 1 


Sij 


= 3cri ■ 


rcTj ■ f 



(2.15) 



l-e-'="'j", (2.16) 

(2.17) 

We note that this Vt^^nr does not contain the part of cTj • cTj-Tj ■ Tj interaction which acquires 
a S{rij) function form in the limit of point particles. This 5— function is probably spread 
out by the finite size of the nucleons, and contributes to the short range part of vnn- How- 
ever, it is difficult to extract it from the phenomenological models. Moreover the dominant 
contribution and the nonlocality effect seem to come from the tensor part of OPEP. 
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Recently, after completion of the present work, Kamada and Glockle [0 found an elegant 
method to obtain a potential VRei that gives exactly the same phase shifts with relativistic 
kinetic energy that a known vnr gives with nonrelativistic kinetic energy. Our objective 
here is not just to find a VRei that is phase equivalent to the Argonne fis; we additionally 
require it to have the Vt^^r^i long range behavior. Both relativistic and nonrelativistic models 
of V contain theoretical long range OPEP; the scattering data is used to determine only the 
phenomenological part vr in these interactions. Our two models are not as exactly phase 
equivalent as Kamada and Glockle's VRei and vnr are, however the differences in their phase 
shifts are negligibly small compared to the uncertainties in the Nijmegen phase shifts ||19|| . 

The parameters of the function Vp{rij) of vr are obtained by fitting the phase shifts and 
deuteron binding energy. Traditionally phase shifts are calculated in configuration space. 



however, in the relativistic case, the Hamiltonian contains i/m^ — which is nonlocal in 
configuration space, therefore we calculate them in momentum space. The details of the 
momentum-space technique have been discussed in Ref. pO[] . 

Some of the important phase shifts are plotted in Fig. |I|. The diamond symbols represent 
the reference nonrelativistic phase shifts obtained with Hj^r, the plus symbols represent 
those calculated from Hr before re-adjusting the parameters in fp(r), and the square symbols 
correspond to those after. The reference phase shifts are almost exactly reproduced by the 
relativistic Hamiltonian Hr as indicated by the good overlap of the diamond and square 
symbols in Fig. |l]. The deviations between the plus and diamond symbols reflect the total 
effect of replacing nonrelativistic kinetic energy and Vj^^nr by the relativistic kinetic energy 
and WjT.ije/ in Hj^r, and are not too large except for the mixing parameter Ei of ^Si —^Di. 
This indicates that relativistic nonlocal effects in two nucleon scattering at E-^^ < 400 MeV 
are rather small. Ei is primarily determined by the tensor force; the relatively large change 
in El is due to the nonlocality of the tensor force in VT^^Rei- 

The new two-body potential is essentially phase equivalent to isoscalar part of Argonne 
f 18 and predicts similar deuteron properties as listed in Table II. Note that the present H^r 
and Hr do not contain electromagnetic interactions. The experimental value of deuteron 
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binding energy (-2.224 MeV) can be obtained from the full Argonne Vig with electromagnetic 
interactions. The 14 operator components of the relativistic and nonrelativistic potentials 
are compared in Fig. ^ Only the static part of the relativistic OPEP, obtained by setting 
m/ \/jP~+lnF equal to unity is used in Fig. |^. Since it is the same as the nonrelativistic 
OPEP, the difference between the potentials shown in the figure is entirely due to that in 
the phenomenological part vr. Equation ( p.5|) shows that the VT^^Rei is smaller than the Vt^^nr 
for p OT p' 0. 

The deuteron S- and D-wave functions are shown in Fig. |^. Here the relativistic D-wave 
is slightly smaller than the nonrelativistic one, presumably because the relativistic tensor 
potential is smaller than the nonrelativistic tensor potential (Eq. |2.5|) for large p and p'. The 
deuteron wave functions in momentum space are shown in Fig. ^. Note that the relativistic 
wave functions are not very different from the nonrelativistic ones. The ratio of the two 
momentum space D-wave functions can be easily understood as discussed below. 

The exact ground state wave function \E' can be expanded in a complete set of states 

\^)=T.'Pi\^)- (2-18) 

i 

For a Hamiltonian given by T + f , where \i) are eigenstates of T, the amplitudes (pi are given 
by 



(^1 


\v\ 


1*) 


(^1 


|T 









(2.19) 



as can be verified from the Schrodinger equation -ff|\E') = -EqI^)- Here Eq is the ground 
state energy, 

In the case of the deuteron we can choose T as the kinetic energy operator and \i) as 
S- and D-waves with momentum p, denoted by \pi), for I = S,D. The amplitudes (f)i{p) 
of these waves give the deuteron wave function in momentum space. We can estimate the 
difference between the nonrelativistic and relativistic deuteron D-state wave function at large 
momentum by assuming that it is primarily generated by the OPEP. In the nonrelativistic 
case this gives: 
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TTl 

(t)D,NR{j)) = :.{pD\Vn,NR\'^), (2.20) 

where we have neglected the Eq in the denominator of Eq. (|2.19|) , since it is much smaller 
than the kinetic energy p^/m at large p. In the relativistic case 



1 TTl 

^D,Rel{p) = / 2 > f T^^=^{PD\Vn,NR\'^), (2.21) 



where the first factor is the relativistic kinetic energy denominator, and the second comes 
from the m/E' factor in the VT^^Rei (Eq. pT3| ). The other m/E factor in the VT^^Rei operates 
on the It is set to unity because most of the deuteron wave function has small relative 
momenta. 

Neglecting the small difference between the relativistic and nonrelativistic the ratio 
of the 4>d{p) is found to be 



^D,Rel{p) P^ 



(2.22) 



4>d,nr{p) 2[^Jm^ + p^ - m)y/m^ + p'^ 

The above estimate is fairly close to the ratio of the calculated D-wave functions as can be 
seen in Fig. ^. Note that this ratio is smaller if the v.,^^fiei is used with the nonrelativistic 
kinetic energy in (pD,Rei{p) (dotted line), and it is larger than one when the Vt,^nr is used 
with the relativistic kinetic energy (dot-dashed line). The relativistic corrections to the 
interaction and kinetic energies have opposite effects on the wave function. The difference 
between the S-wave functions is influenced by the changes in the kinetic energy and the vr. 
The effects of these changes on the phase shifts and the deuteron energy must cancel by 
construction, and they seem to largely cancel in the (ps{p)- 

Some of the deuteron momentum space results are listed in Table III. This table offers a 
microscopic view of how various relativistic effects were buried in the nonrelativistic models. 
Relativistic nonlocalities reduce the OPEP contribution by ~2.6 MeV, while the relativistic 
kinetic energy is smaller by ~ — 1 MeV, giving a net effect of 1.6 MeV which is canceled by 
the change in the phenomenological vji. 

The variational Monte Carlo calculations for the A > 2 systems have to be carried out 
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in configuration space. We therefore liave to Fourier transform the Vt^^rsI which depends on 
both p and p' (Eq. p. 51), yielding a nonlocal potential in configuration space: 



v.Mr\ r) = / e--'-' p) (2.23) 

The exact integral in Eq. (|2.23| ) is extremely difficult to calculate. The series obtained by 
expanding Vn,Rei{p', p) in powers of /m^ is given by: 

However, this series does not have good convergence. In the case of deuteron, the expectation 
value of Vjr^NR, for the eigenfunction of our relativistic Hamiltonian is -21.39 MeV. The term 
in Vir^Rei, of order 1/m^, contributes 3.48 MeV to the {vTT,Rei), while that of order l/m"^ gives 
-1.45 MeV, and the exact {v^^Rei — Vtt^nr) is 2.59 MeV. Therefore the series converges slowly 
to the exact value. This may appear surprising because the expectation value of the kinetic 
energy of the deuteron (Table III) is only about 20 MeV, giving /m? k, 0.02 on average. 
However, the deuteron has large momentum components via the D-wave, or equivalently 
the tensor correlations, and most of the OPEP contribution is from those. Thus it is not 
surprising that the expectation value of OPEP is sensitive to higher powers of nucleon 
velocities. 

For the relativistic OPEP, a good convergence is achieved by using the variables: 

Q = ^(p + p'), q = p-p', (2.25) 

x=^(r + r'), y = r-r', (2.26) 

for which 

2 

Tfl 

Vn,Rel = Vn,NR{q) 



Q' + i ,1 m2(Q.q) 



1- „ + 



Here we expanded the Wtt.m in powers of (Q • q)^/(m^ + Q'^ + g^/4)^. This series appears 
to converge rapidly. In the case of deuteron, the leading relativistic correction given by 



(2.27) 
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the second term is 2.7 MeV, the third term gives -0.18 MeV, and the exact value is 2.59 
MeV. The third term contains Oqq dependence and results in complicated operator forms as 
shown in the Appendix. Moreover, the third and higher terms account for only ~ 4 % of the 
relativistic correction to OPEP expectation value in the deuteron {i.e. ~ 0.6 % of (v-^^Rei))- 
Therefore only the first two terms are considered in this work, and the relativistic OPEP is 
approximated by: 

,2 



2 ^^TT.Affliqj- 



(2.28) 



With this w,r,_Re« in Eq. ( p.9|) , we refit the phase shifts and deuteron binding energy. The 
parameters in Vp{r) are very similar to those obtained with the exact v^^^Rei given by Eq. ( p. 51) . 
The new phase shifts and the relativistic potentials are very similar to those shown in Figs. 
|1] and |. 

The configuration space potential (Eq. |2.23|) is given by: 



W7r,MlX,y) 



d^Q d?q 



m 



,«(Qy+q-x) 



(27r)3 (27r)3 ^2 + g2 + ^ 
and is simple to evaluate. The v.„^]^r (Eq. |2.14| ) in momentum space is given by: 



^ /2 



y^{q) (Ti ■ (Tj + %{q) (Ti ■ (Tj -(Ti ■ qcTj • q 



where 



(2.29) 



(2.30) 



yn{q) = J Y^{r)jo{qr)r^dr, 
%{q) = J T^{r)j2{qr)r'^dr. 



Substituting these into Eq. ( p.29| ) gives 



Vn,Rel{x, y) = 3/^^:^ fiv) [F<tt{x, y) (Ti ■ (T j + Ftr{x, y) Sij{x, X)] Ti ■ Tj, 



(2.31) 
(2.32) 



(2.33) 



with 
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F.r{x, y) = - I q'dq y^q) Uqx) e-(V-^+''V4-n^) ^2.34) 

FUx, y) = - j q'dq %{q) Uqx) e-(V-^+«V4-™) ^2.35) 
, rn? e~^y 

fiy) = J , 2.36 

Sij{x, X) = 3 CTi • X CTj • X — CTj ■ (Tj. (2.37) 

In the limit m oo, f{y) becomes 5^(y), F„r{x,y) Yt,{x) and Ftr^x.y) T^(a;). When 
y — > 0, we have r' = r, x = r and Eq. (|2.33|) becomes Vt^^nr. Figure | shows F„r{x,y) 
and Ftr{x,y) as a function of x for various values of y. Note that the solid lines for y = 
correspond to the nonrelativistic Yt^{x) and T^(x). The volume integral of F„T-{x,y) is 
independent of y, whereas that of Ftr{x,y) decreases with y. Therefore relativistic effects 
mostly come from the tensor part of OPEP. 

III. VARIATIONAL MONTE CARLO CALCULATIONS AND RESULTS 

A. VMC techniques 

With the relativistic Hamiltonian discussed in the previous section, we can proceed to 
evaluate the energy expectation value 

(Hr) = (T) + (v,,) + {6v,,) + (3.1) 

for A > 3 nuclei using the Monte Carlo technique. 

The Monte Carlo method [^] offers a useful way to handle the multidimensional integrals 
which would otherwise be impractical by the usual numerical methods. The basis of this 
method is that instead of integrating over a regular array of points, we sum over a set of 
configurations {Ri} distributed with probability w(R). Here R = (ri, r2, ■ ■ ■ , r^) denotes 
the configuration of all the nucleons in the nucleus. There are various techniques for sampling 
w(R) and in this work Metropolis sampling method is used to treat the complicated 
distributions. 
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Variational Monte Carlo (VMC) technique is based on variational principle that the 
minimum expectation value of the Hamiltonian is closest to the ground state energy of the 
system. Starting from a variational wave function, which depends upon several variational 
parameters {ai, q;2, • • • , we evaluate the expectation value of the Hamiltonian using the 
Monte Carlo configuration samples R,: 

/dR^t(R)^^„(R) 



/dR^l(R)^„(R) 

±d, (3.2) 



where S is the standard deviation. Typically block averaging scheme is used to obtain a 
normal distribution and the error can be conveniently evaluated from it. We divide Nc 
configurations into Ni, blocks each containing A^o = ^c/^b configurations. The average 

is evaluated for each block. The expectation value of H is given by 

1 ^6 

W^^Y^Hb (3.4) 

with the standard deviation 



(3.6) 



The Monte Carlo result is exact when the number of configurations Nc ^ oo, although in 
practice A^c=50000 seems to be enough to obtain results with sufficiently small statistical 
errors. The weight function is usually chosen to be 

w{Ri) = *|;(R,)*.(RO (3.6) 

to maintain small Monte Carlo error. Note that when "i/y is the eigenstate of H, the Monte 
Carlo sampling error becomes zero. Finally, the parameters (cti, a2, - ■ ■ , an) are varied to 
minimize the energy. 
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Some of the terms in (Hr) (Eq. can be calculated straightforwardly and have been 
discussed in Refs. and ||T^ . The terms that require special techniques are the relativistic 
kinetic energy {^/m"^ — V^) and VT^^Rei in the two-body potential. The kinetic energy term 



has been calculated previously in Ref. and the calculation of (v-^^Rei) is discussed in Sec. 

inc. 

B. Relativistic wave functions 

In the nonrelativistic case, variational wave functions of the form 

1^.)= (l+ E Fv>] lsl[F^^ m, (3.7) 

\ i<j<k J \ i<j J 

having symmetrized product of pair correlation operators Fij and a sum of triplet correlations 
Fijk operating on an antisymmetric, uncorrelated wave function |$), have been commonly 
used. The Fij and Fijk correlation operators reflect the effects of two-, three-body interac- 
tions on the wave function. The uncorrelated wave function has no spatial dependence for 
A <4 nuclei. A good representation of such a wave function is given in Ref. . 

The pair correlation operator F^j is constructed from correlation functions which satisfy 
Schrodinger-like 2-body equations, with appropriate boundary conditions. Their solutions 
are like deuteron wave functions "^nr and displayed in Fig. In the case of A=2 
deuteron, both nonrelativistic and relativistic correlation functions can be easily solved in 
momentum space; they are not very different from each other as can be seen in Fig. ^. In 
A >2 nuclei, the nonrelativistic pair correlation equations can be easily solved in configu- 
rations space, however, the relativistic equations are more difficult to solve. Therefore we 
seek good approximations for the relativistic pair correlation functions. 

Our method can be easily illustrated using the example of deuteron. Its variational wave 
function is expressed as: 

^ = ^NR + K'^r-'^nr), (3.8) 
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is calculated using VMC, and A is varied to minimize it. The results are shown in 
Fig. 1^. The error bars shown in Fig. ^ originate from the statistical sampling and are <1% 
of the binding energy. The same configurations are used to calculate the energies for all A, 
hence the errors are correlated. The minimum value of (Hfi) does occur at A = 1 where 
\I' = "^fi as expected. The difference in (Hji) between A = (using a nonrelativistic wave 
function) and A = 1 (using a relativistic one) is ~0.04 MeV. This means that if we were 
to use the nonrelativistic wave function to calculate the expectation value of Hn, the result 
will be off by only 2% for the deuteron. 

In heavier nuclei we also expect the optimum nonrelativistic wave function to provide 
a good approximation for relativistic wave function. The difference between the two is 



define: 





and ^Sq 


correlation functions at small r. 


We therefore 


/o,i — 


fo,l,NR ( 




(3.9) 


fi,o ~ 


flfi,NR ( 


;i + A^5j, 


(3.10) 


/l,0 ~ 


fto,NR ( 


;i + a^dJ, 


(3.11) 



where ior the channel c is defined as: 

Ur) = ^-^^'^ - (3.12) 

<Pc,NR{r) 

At small r, the central /^q and tensor /^q correlation functions in spin-isospin S,T = 1,0 
states are proportional to the S and D radial wave functions of the deuteron, respectively. 
Therefore the ^35^ and C,iDi can be calculated exactly in momentum space using the rela- 
tivistic and nonrelativistic deuteron wave functions for the 0c,_r and (p^NR- These ^'s are 
rather short-ranged (Fig. |^), and we do not expect them to vary significantly in larger nuclei. 
Wiringa has shown that the nonrelativistic correlation functions for ^H, and ^He are 
almost the same at small r. 

A similar calculation of ^i^^ is not possible because of the absence of a bound state in 
that channel. However, we can obtain an artificial bound state by slightly increasing 
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the strength of intermediate range attraction in v{^So). The binding energy and the wave 
function at large r are very sensitive to small changes in v{^Sq), however, ,^15^^ is relatively 
insensitive. As an example, the ^i^^ obtained from artificial ^Sq bound states with energies 
of -1 and -10 MeV, shown in Fig. |^, are very similar. 

The VMC energy of with the relativistic Hamiltonian is shown as a function of A 
in Fig. p. The minimum occurs at A=0.5 instead of 1 (the expected value), however, the 
difference in (Hn) between A=0.5 and 1 is rather small and of the order of the Monte Carlo 
sampling error. The minimum energy for ^He occurs at the expected A=1.0. 



C. Expectation value of the nonlocal potential 

Consider an A-nucleon system whose wave function is denoted as "^{ri, r2, ■ ■ ■ , r^). We 
define: 

x, = ^(r, + r:), x, = l(r,. + r;), (3.13) 

so that 



x = Xi-Xj, ^ = ri-r'i = r'j-rj, (3.14) 



as illustrated in Fig. 10. The expectation value of relativistic OPEP is then given by 



{vn,Rei) = E / n d^Tkd^Xid^Xjd^y f ri, ■ ■ ■ , Xi - ^, ■ ■ ■ , Xj- + ^, ■ ■ ■ , J 

i<j'' k^i,j ^ 4 4/ 

^u,fe«(|xi-Xj|,?/) ^ (^ri,---,Xi + |,---,Xj - |,---,r^^ , (3.15) 

where VTr,Rei (|xi — Xj|,?/) is previously calculated in Eq. ( |2.33| ). The integration over the 
Tfc's, Xj, Xj and the solid angle of y is carried out by the Monte Carlo method, while that 
over the magnitude of y is carried out with Gauss-Laguerre integral. 



D. VMC results 



The VMC results for '^H and '^He are listed in Table IV. Note that in principle we should 
use GFMC to calculate the exact binding energies, but the relativistic effects resulting from 
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the difference between {Hji) and (Hnr) are small and presumably not too different from 
those estimated using VMC. 

The total relativistic effect on the binding energy is ~0.3 MeV for and ~1.8 MeV 
for ^He. Most of the effect comes from the boost correction which is 0.42 MeV for '^H and 
1.94 MeV for ^He. The net effect of relativistic corrections to the kinetic energy and the 
two-body potential, on the binding energy is rather small: ~ —0.12 ±0.06 MeV in (^H) and 
~ —0.17 ± 0.10 MeV in (^He). Since both Hn^ and are constrained to give the same 
deuteron binding energy, the changes in (T) and {vij) cancel exactly in ^H. In and ^He 
they appear to largely cancel and give a rather small net effect. 

In view of the slow convergence of the {IvT^^Rei — "WTr.Afj?!), when expanded in powers of 
/rn?, one may question the validity of calculating the boost interaction 5v only up to 
first order in P'^/Am?. The average kinetic energies of nucleons in nuclei are rather small 
giving average pj/m? < 0.1. The expansion has good convergence for such values. However, 
two nucleons can occasionally have large relative momenta when they come close together 
as illustrated in Fig. In that configuration ~ p and Pj ~ — p are both large due to 
strong short-range interaction between nucleons i and j. Such configurations are responsible 
for the slow convergence of the expansion of {\vT,^Rei — Vt,^nr\)- In these configurations the 
squares of the total pair momenta are of order Pj| ~ Pjj. ~ p^, while Pf^ has small, near 
average value. Thus the expansion parameter P^/Am^ for the boost interaction is effectively 
four times smaller than that of {v^^^Rei ~ Vj^^nr) when the momenta are generated by pair 
correlations, therefore we expect the boost expansion to converge more rapidly. Moreover, 
{\vTT,Rei — Vt,^nr\) is much larger (~ 6 and 13 MeV in and ''He respectively) than {\5v\) 
and has to be calculated with higher relative accuracy to obtain a final total energy with 
error ~ 1 %. 
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IV. CONCLUSIONS AND OUTLOOK 



We find that the relativistic effects in on-shell OPEP are quite substantial. The expecta- 
tion values of 'y,r,_Re« are smaller than those of Vn,NR by ~ 15 %. Since the expectation values 
of OPEP are much larger than nuclear binding energies, the differences in the OPEP expec- 
tation values are comparable to the total nuclear energy. However, nuclear Hamiltonians are 
not derived from first principles, they are obtained by fitting data. The substantial difference 
between v-^^nei and v-„^nr is compensated in the Hnji by that in the kinetic energy T and vr 
so that it gives the same scattering cross sections and deuteron energy as the Hr. We find 
that this compensation works rather well for three- and four-body nuclei. In absence of the 
boost interaction our nonrelativistic and relativistic Hamiltonians seem to give very similar 
results for the binding energies and wave functions of light nuclei. It is probably necessary 
to examine one-pion exchange current contributions to elastic scattering form factors and 
radiative capture reactions to see the effect of the m/E factors in the OPEP. 

The modern two-nucleon potential models, which include the Nijmegen models I, II and 
Reid-93 [Q, Argonne vis and CD-Bonn accurately reproduce the NN-scattering 



data in the Nijmegen data base. Friar et al. [|28| have studied the triton energy with the 
Nijmegen and Argonne models, without boost or three-nucleon interactions, using accurate 
Faddeev calculations. The energies obtained with the three local potential models, Reid-93, 
Nijmegen-II and Argonne fig are respectively -7.63, -7.62 and -7.61 MeV. These energies 
are very close, and these models also give very similar values (5.70, 5.64 and 5.76 %) for 
Pdi the fraction of D-state in the deuteron. The boson exchange Nijmegen-I model contains 
nonlocal terms and gives -7.72 MeV for triton energy and 5.66% for P^. Comparison of the 
results of Nijmegen I and II models indicates that the total effect of the nonlocalities on 
energies and the wave functions could be small. The present results support this conclusion; 
inclusion of relativistic nonlocalities in on-shell OPEP and kinetic energy lowers the triton 
energy by ~ 0.1 MeV and Pd by 0.04%. 

In contrast the CD-Bonn potential gives rather different results from the Nijmegen and 
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Argonne models. It gives a triton energy of -8.00 MeV and Pd = 4.83%. The OPEP in 
the CD-Bonn model has additional off-shell nonlocalities which disappear in the on-shell 



relativistic OPEP. It is defined as 29 



V7T,CDB[P,P) 



fl 



CTj ■ q cr ■ q 



-{E' - E) 



'(Ti • p CTj • p CTj • p' (Tj ■ p'' 



(4.1) 



E + m E' + m 

where E = ■\/m^~+~p^, E' = y/ m'^ + p''^ . The term proportional to {E' — E) does not 
contribute to the on-shell OPEP, and is absent from our "VT^^Rei given by Eq. ( p.5| ). The 
OPEP gives large contributions by coupling states with small p to large p'. Therefore we 
consider the case p = for which p' = — q, and expand in powers of /rn?. This gives: 



,NR 



1 - 



3g2 
4m2 

2m2 



+ ■ 



+ ■ 



(4.2) 
(4.3) 



indicating presence of larger relativistic corrections in the Vt^^cdb- 

The second order contributions of the Vt^^cdb provide a sightly better approximation 
to the sum of the twelve time-ordered two-pion exchange diagrams p2| with only positive 
energy nucleons in the intermediate states, in relativistic field theories with pseudoscalar 
coupling. The two-pion exchange diagrams with antinucleons are discarded on arguments 
based on chiral symmetry. It is not obvious that this off-shell term must be retained, and 
not discarded. It is necessary to find experimental tests for its existence, as well as for the 
suppression of OPEP hj m/E factors considered in this work. 

The off-shell behavior of the OPEP can also be changed by using combinations of pseu- 
doscalar and pseudovector couplings. In Friar's notation the possible off-shell behaviors 
are characterized with parameters /2 and z/, and up to order p'^/m^ they are related by uni- 
tary transformations. Up to this order our v-jt^Rei has an off-shell behavior with u = 1/2 
and /i = 0, while that of v-^^cdb has v = 1/2 and /i = —1. The Pd has smaller values for 



fi = —1 used in the CD-Bonn potential. 
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The boost interaction Sv gives the dominant relativistic correction to the binding energies 
of hght nuclei in the present formahsm. Contributions of 6v are very small, only ~ 1 % of 
that of V in and ^He. However they are not canceled by relativistic effects in either T or v, 
and therefore dominate the net effect. Like the two-nucleon, the three-nucleon interaction 
is also not derived from first principles. Urbana models of Vijk contain two terms: the 
attractive two-pion exchange term V-^l and a repulsive phenomenological term V^jf.. Their 
strengths are chosen to reproduce the triton energy and the density of nuclear matter without 
considering any relativistic effects. In light nuclei the repulsive 6v contribution is about 37 
% of that of Vj^j.. Thus the strength of V^^f. in Hr has to be reduced by 37 % to obtain 
the experimental energies of light nuclei. One could also choose to increase the strength of 
V^^l or some combination of increasing V^^l and decreasing V^^^, but in this work we keep 
V^^l unchanged. The difference in the three nucleon interaction then compensates for the 
omission of 6v in conventional nonrelativistic nuclear Hamiltonians. It appears that this 
compensation works rather well in light nuclei having up to eight nucleons as well as in 



nuclear and neutron matter up to normal densities |^ . It presumably works also in heavier 
nuclei where it is not tested. However, at several times nuclear matter densities, encountered 
in neutron stars, the effective three-nucleon interaction overestimates the 6v contribution 
significantly 13^ . 
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APPENDIX A: RELATIVISTIC OPEP IN CONFIGURATION SPACE 



The third term in Eq. ( 2.27 ) gives the second-order relativistic correction to OPEP and 
is denoted as v^"^^: 



^(2)(x,y) 



(27r)3 (27r)3 

where 9qQ is the angle between q and Q, and 

1 m?Q'^q^ 



Expressing cos^ O^q as 



9{<1,Q) 



2 (m2 + g2 + 



2l^3 

4 



COS^^qQ = -P2(cOSe5Q) + -, 



(Al) 



(A2) 



(A3) 



we get 



+ 



3 J (27r)3(27r)3 
2 f d^Q d^q 



g{q, Q) P2(cose,Q) v^,nr{^) e^(Q 



(A4) 



3 J (27r)3(27r)3 

(2) 

The first integral denoted by ^1 (x, y) is independent of 9qQ and can be easily evaluated by 
using Eq. ( p.SOj ) and the following identities: 



Im 

Yi*^{q)Yi,m'{q)dVtq = Su'Smm', 

v{r)(Ti ■ f(Tj ■ r e^^'^d^r = 



(Ti ■ Vqcr 



We obtain: 



(2)/ N 

^1 (x,y) 



^^ fnNN 



Q^dQq'^dqg{q,Q)y^{q) jo{Qy)jo{qx) (Ti ■ crj 



+ Q dQq dq g{q,Q) %{q) jo{Qy)j2{qx) Sij{x, x) 



Ti ■ Ti 



where 3^7r(g) and %{q) are previously given in Eqs. ( ^^.311 ) and ( |2.3^ 
To calculate the second integral in Eq. we use 



(A5) 

(A6) 
(A7) 



(A8) 
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A'TT 

The integral over the sohd angles becomes: 

dfiQdfi,e*(Q-y+'i-'')P2(cos^,Q) = {A7rfP2{cose,y)MQy)j2{qx), 
and the second term in Eq. (|A^), denoted by 1^2^'' (x, y) is: 



(A9) 



(AlO) 



97r3 An 



JP2(cos^,,) jQMQqMqg{q,Q) [y^q) + Mq)] J2{Qy)j2iqx)cT,-(Tj, 
P2{cose,y)jQMQdqg{q,Q)%{q)UQy)j2{qx) }T,-r,. (All) 



Here O^y is the angle between x and y. The gradient operators in the second term act 
on both P2(cos6'^.j,) and j2{qx). 

The Q-integral in Eqs. ( |A8| ) and ( [All| ) can be performed analytically and results 



J Q'dQg{q,Q)UQy) = f{y) [Z,{q,y) - Z^{q,y))] 
Q'dQg{q,Q)UQy) = fiy)Z^iq,y), 



where f(y) is the Yukawa function given in Eq.(2.33) and 



Z2{q,y) 



m? + ^ y 



Ziiq,y). 



v^"^^ is finally obtained as: 



(2)^„ _ fnNN 



3 47r 



f{y)[Ii{x,y,9„j) (Ti ■ (Tj + l2{x,y,9^y) Sij{x,x) + 

hix, y, 6^y) Sijiy, y) + h{x, y, O^ry) Sij{x, y)]Ti ■ tj, 



(A12) 
(A13) 



(AM) 
(A15) 



(A16) 



where the tensor operator Sij{x,y) is defined as: 



Sij{x, y) = 2 i^i ■ ^^j ■y + o-i-y(Tj-x)-x-ycri-(Tj, 



and Ji, I2, h and I4 are given by: 



(A17) 
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h{x, y, 9^y) = J^^oix, y) - J^f^ix, y) + 2 J^f^i^, y) P2(cos 

hix, y, e^y) = J^22{x, y) - 3 Hi^iix, y) + 'iJ^l^ix, y) cos^ 
h{x,y,6^y) = &nl2{x,y) 
h{x,y,9xy) = -I2nls{x,y) cosO^y 

nli{x,y) 



1 



37r3 

1 



q'^dqZa{q,y)yn{q) Ji{qx) 

3^3 y q^(iqZa(q,y)%(q) ji(qx) 
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FIGURES 

FIG. 1. Phase shifts as a function of lab energy, o: the reference phase shifts obtained with 
Hnr] +'■ those using vr from H^n in H^; □: those with the new relativistic Hamiltonian Hn 
with re-adjusted vr. 

FIG. 2. Comparison of relativistic (solid lines) and nonrelativistic (dashed lines) potentials 
vi — f 14 of the operator format (refer to Eq. |2.12D . Note that f 4 and f g contain contributions from 
both and vr, and only a local nonrelativistic OPEP is used. 

FIG. 3. Deuteron wave functions. 

FIG. 4. Deuteron wave functions in momentum space. 

FIG. 5. Ratio of deuteron relativistic to nonrelativistic D-wave function in momentum space 
(solid line), and the simple estimate of Eq. ^.22| (dashed line). The dotted and dot-dashed lines 
represent results calculated using Tnr-, Vjr^Rei and Trsi, v-n^NR in 4>D,Rei in Eq. ( 2.22| ), respectively. 



The wiggle in the calculated ratio around 7 fm comes from a node in the wave function. 



FIG. 6. Fg-r and as a function of x for various values of y. 



FIG. 7. VMC results for deuteron with 100,000 configurations. 



FIG. 8. ^ = ^^jj^ as a function of r. 



FIG. 9. VMC results for triton with 50,000 configurations. 



FIG. 10. A diagram to illustrate the calculation of nonlocal interaction contributions. 
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FIG. 11. A naive picture of to illustrate large momentum contribution from a configuration 
where two nucleons are close together. 
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TABLES 

TABLE I. Nonrelativistic Green's Function Monte Carlo (GFMC) results (in MeV) for light 
nuclei with the Argonne vis and Urbana IX potentials. The first line gives the experimental energy 
while the next four list the calculated total, kinetic, two- and three-body interaction energies. The 
last two lines give the contribution of the pion exchange parts of 2- and 3-body interactions. 





2H 


3h 


%e 


6Li 


^Li 


^Be 


E 


-2.2246 


-8.48 


-28.30 


-31.99 


-39.24 


-56.50 


<E> 


-2.2248(5) 


-8.47(1) 


-28.30(2) 


-31.25(11) 


-37.44(28) 


-54.66(64) 


<T> 


19.81 


50.8(5) 


111.9(6) 


150.8(10) 


186.4(28) 


246.3(56) 


< Vij > 


-22.04 


-58.4(5) 


-135.4(6) 


-179.2(10) 


-220.8(30) 


-295.8(62) 


< Vijk > 


0.0 


-1.20(2) 


-6.4(1) 


-7.2(1) 


-8.9(2) 


-14.8(5) 




-21.28 


-43.8(2) 


-99.4(2) 


-128.9(5) 


-152.5(7) 


-224.1(9) 


< V^^ > 


0.0 


-2.17(1) 


-11.7(1) 


-13.5(1) 


-17.1(4) 


-28.1(8) 



TABLE IL Deuteron properties 

Hnr Hr 

binding energy (McV)'^ -2.242 -2.242 

quadrupole moment (fm^) 0.269 0.271 

% of D state (Pd) 5.776 5.732 
t without electromagnetic interactions 
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TABLE III. Results of momentum space deuteron calculations. 



{^nr\Hnr\^nr) {^r\Hr\-^r) 



{E) -2.242 -2.242 

(T) 19.882 18.877 

{vij) -22.125 -21.119 

{v^) -21.356 -18.797 

{vr) = {vij - v„) -0.769 -2.322 

{Vn,Rel - Vn,NR) 2.589 



TABLE IV. VMC results for and ^He, calculated with 50,000 configurations. 
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